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I.  INTRODUCTION 


In  the  summer  of  1974,  a  155mm  shell  with  a  liquid  payload  was  fired  with 
a  yawsonde  telemetry  unit  and  this  projectile  developed  a  large  amplitude  con¬ 
ing  motion.  When  the  coning  motion  exceeded  forty  degrees,  the  telemetry 
record  showed  a  very  rapid  despin  of  forty  percent  in  less  than  five  seconds.1 
Since  1974,  ten  more  observations  of  liquid-filled  projectiles  with  large  con¬ 
ing  motions  and  large  despin  rates  have  been  made.2  4  In  addition,  thirty- 

six  observations  have  been  made  of  large  coning  motions  and  despin  moments  for 
projectiles  with  payloads  of  liquid-solid  mixtures.5  Indeed,  all  observed 
flights  of  projectiles  with  liquid  or  liquid-solid  payloads  that  performed 
large  coning  motion  also  experienced  large  losses  in  spin. 

This  characteri Stic  association  of  a  large  despin  moment  with  a  large 
coning  motion  for  a  projectile  with  a  moving  payload  can  become  a  diagnostic 
tool.  Miller6  suggested  using  this  tool  to  determine  the  existence  of  small- 
amplitude  unstable  liquid-induced  side  moments.  Spin  measurements  made  during 
coning  motion  on  a  gyroscope  predicted  flight  pitch  instabilities  caused  by 
very  viscous  liquids,  and  these  were  observed  in  flight.7 


1.  D'Amico,  W.  P. ,  Oskay,  V .,  and  Clay ,  W.  H.,  "Flight  Tests  of  the  155mm 
XM687  Binary  Projectile  and  Associated  Design  Modification  Prior  to  the 
Nicolet  Winter  Test  1974-197 5 US  Ballistic  Research  Laboratory ,  Aberdeen 

Proving  Ground ,  Maryland ,  BRL  Memorandum  Report  ARBRL-MR- 2748,  May  1977. 
(AD  B019969). 

2.  D'Amico,  W.  P.,  Clay ,  W .  H.,  and  Mark,  A.,  "Yawsonde  Data  for  M687-Type 
Projectiles  with  Application  to  Rapid  Spin  Decay  and  Stewartson-Type  Spin- 
Up  Instabilities ,"  US  Ballistic  Research  Laboratory ,  Aberdeen  Proving 
Ground,  Maryland,  BRL  Memorandum  Report  ARBRL-MR-03027 ,  June  1980. 

(AD  AO 89 6 46) 

3.  D'Amico,  W.  P.,  Clay ,  W.  H.,  and  Mark,  A.,  "Diagnostic  Tests  for  Wick-Type 
Payloads  and  High-Viscosity  Liquids,"  US  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland,  BRL  Memorandum  Report  ARBRL-MR-2913 , 
April  1979.  (AD  A072812). 

4.  D'Amico,  W.  P. ,  and  Clay,  W.  H.,  "High  Viscosity  Liquid  Payload  Yawsonde 
Data  for  Small  Launch  Yaws,"  US  Ballistic  Research  Laboratory ,  Aberdeen 

Proving  Ground,  Maryland,  BRL  Memorandum  Report  ARBRL-03029 ,  June  1980. 

(AD  AO 8 84 11) . 

5.  Murphy,  C.  H.,  "Liquid  Payload  Roll  Moment  Induced  by  a  Spinning  and 
Coning  Projectile ,"  US  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground,  Maryland,  BRL  Technical  Report  ARBRL-02521 ,  September  1983  (AD 

A133681).  (See  also  AIAA-83-2142 ,  August  1983.) 

6.  Miller,  M.  C .,  "Flight  Instabilities  of  Spinning  Projectiles  Having 
Nonrigid  Payloads ,"  Journal  of  Guidance,  Control,  and  Dynamics,  Vol.  5, 
March-April  1982,  pp.  151-157. 

7.  D'Amico ,  W.  P. ,  and  Miller,  M.  C.,  "Flight  Instability  Produced  by  a 
Rapidly  Spinning ,  Highly  Viscous  Liquid ,"  Journal  of  Spacecraft  and 
Rockets,  Vol.  16,  January -February  1979,  pp.  62-64. 


5 


The  linear  liquid-induced  side  moment  was  first  computed  by  Stewartson8 
for  an  inviscid  liquid  payload  by  use  of  eigenfrequencies  determined  by  the 
fineness  ratio  of  the  cylindrical  container.  Wedemeyer9  introduced  boundary 
layers  on  the  walls  of  the  container  and  was  able  to  determine  viscous  correc¬ 
tions  for  Stewartson's  eigenfrequencies,  which  could  then  be  used  in  Stewart- 
son's  side  moment  calculation.  Murphy10  then  completed  the  linear  boundary 
layer  theory  by  including  all  pressure  and  wall  shear  contributions  to  the 
liquid-induced  side  moment. 

The  first  theoretical  work  on  liquid-induced  roll  moments  was  done  by 
Vaughn11  in  1978.  Although  fair  agreement  with  Miller's  data  was  obtained, 
the  work  was  marred  by  some  hard-to-justi fy  algebraic  steps.  Recently  Vaughn 
et  al  12  developed  an  impressive  numerical  capability  to  compute  roll  moments 
for  very  high  viscosity  liquids  and  obtained  excellent  agreement  with  Miller's 
data.  In  Reference  5  the  linearized  Navier-Stokes  equations  were  used  to 
develop  a  relationship  between  the  liquid  side  moment  and  the  liquid  roll 
moment.  It  is  the  purpose  of  this  note  to  give  the  results  of  this  analysis 
and  illustrate  its  predictions  with  yawsonde  data  recently  analyzed  by  Pope.13 


II.  LIQUID  ROLL  MOMENT 

The  pitching  and  yawing  motion  of  a  symmetric  spinning  projectile  can  be 
represented  as  the  sum  of  two  rotating  exponentially  growing  two-dimensional 
vectors.  In  terms  of  complex  variables,  this  relationship  assumes  the  form 

^  i  <h  i  4>o 

5  =  3  +  i  a  =  Kje  +  K2e  £  (2.1) 


8.  Stewartson ,  K. ,  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid," 
Journal  of  Fluid  Mechanics,  Vol.  5 ,  Part  4 ,  September  1959 ,  pp.  577-592. 

9.  Wedemeyer ,  E.  H. ,  "Viscous  Correction  to  Stewartson's  Stability 
Criterion ,"  US  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground , 
Maryland ,  BRL  Report  1325,  June  1966.  (AD  489687). 

10.  Murphy,  C.  H .,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous 
Liquid  Payload ,"  US  Ballistic  Research  Laboratory ,  Aberdeen  Proving 
Ground ,  Maryland,  BRL  Memorandum  Report  ARBRL-MR-03194 ,  August  1982. 

(AD  A118676) .  (See  also  Journal  of  Guidance,  Control,  and  Dynamics,  Vol. 
6,  July-August  1983 ,  pp.  280^2867)  " 

11.  Vaughn,  H.  R.,  "Flight  Dynamic  Instabilities  of  Fluid-Filled  Projec¬ 
tiles,"  Sandia  Laboratories ,  Albuquerque,  NM,  SAND  78-0999 ,  June  1978. 

12.  Vaughn,  H.  R. ,  Oberkampf ,  W.  L. ,  and  Wolfe,  W.  P. ,  "Numerical  Solution 
for  a  Spinning,  Nutating,  Fluid-Filled  Cylinder ,"  Sandia  Laboratories , 
Albuquerque,  NM,  SAND  83-1789,  December  1983. 

13.  Pope,  R.  L. ,  "Further  Analysis  of  Yawsonde  Data  from  Some  Liquid  Payload 
Projectiles ,"  US  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground, 
Maryland,  BRL  Memorandum  Report  03329,  December  1983.  (AD  A137256 ) 
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where  In  (Kj/Kjg)  =  Gj  tj  <Jt 
+j  =  ^jO  +  Tj^ 

$  is  the  roll  rate  and  is  assumed  to  be  positive. 


Two  coordinate  systems  will  be  used  in  this  note:  the  nonrolling  aero- 

ballistic  XYZ  system  whose  X  -axis  is  fixed  along  the  missile's  axis  of 
symmetry  and  the  inertial  XYZ  system  whose  X-axis  is  tangent  to  the  trajectory 
at  time  zero.  Both  coordinate  systems  have  origins  at  the  center  of  the 
cylindrical  payload  cavity,  which  is  assumed  to  be  at  the  center  of  mass  of 
the  projectile. 


The  transverse  moment  exerted  by  a  liquid  payload  was  expressed  in 
Reference  10  in  the  following  form. 


MLY  +  1  MLZ  "  mL  a2^2  ^-Tl  CLM1  Kle  +  x2  CLM2  K2e  ^ 


(2.2) 


where 


ni|_  =  mass  of  liquid  in  fully  filled  container 


and 


a  =  maximum  radius  of  liquid  container. 


For  linear  fluid  motion, 
fill  ratio,  and  shape  of 


CLMl  should  depend  on  , 
cavity.  A  similar  remark 


>  time,  Reynolds  number, 
applies  to  C|_^. 


The  C|_M.  are  complex  quantities  whose  imaginary  parts  represent  in-plane 
moments  causing  rotation  in  the  plane  of  exp  (i<f>.)  and  whose  real  parts  repre- 

J 

sent  side  moments  causing  rotation  out  of  the  plane  of  exp  (i <j>-)»  Thus,  C|_^. 
can  be  written  as  J  J 

CLM .  =  CLSM .  +  1  CLIM .  (2*3) 

J  J  J 

where  Cj_g^ .  are  liquid  side  moment  coefficients  and  are  liquid  in-plane 

J  J 

moment  coefficients. 


The  liquid  roll  moment  can  now  be  defined  in  a  similar  way.  Symmetry 
considerations  imply  an  even  dependence  on  coning  amplitudes. 
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(2.4) 


MLX 


a2*2 


[C 


LRMr 


rl  K1  CLRM,  +  t2  K2  CLRM0^ 


where  are  functions  of  Reynolds  number,  fill  ratio,  shape  of  cavity,  and 

J 

time.  In  addition,  for  j  =  1  and  2,  the  Clr^.  may  be  functions  of  t^,  ej  and 

2  •  ^ 

Kj  .  is  the  roll  moment  coefficient  exerted  by  the  liquid  during  spin- 

up  and  is  zero  when  steady  state  is  reached.  For  the  steady-state  coning 
motion  considered  in  the  Stewartson-Wedemeyer  theory  (C^r^  =<2=0),  only 


one  roll  moment  coefficient  is  present  and  we  will  omit  the  subscript  on 


and  C 


LRM-, 


m  -  -  m  a2l2  v  2  r 
Mlx  mL  a  <f>  x  K1  CLRM. 


(2.5) 


Throughout  this  note,  we  will  consider  only  roll  moments  that  can  be  approxi¬ 
mated  by  the  one-term  Eq.  (2.5)  rather  than  the  more  complete  three-term  Eq. 
(2.4).  Similarly,  the  transverse  liquid  moment  will  be  approximated  by  the 
first  term  of  Eq.  (2.2). 


M|^y  ^  i  M^^ 


a2-2 
a  <(>  t 


(^LSM  +  tm) 


LIM' 


Kxe 


<t>l 


(2.6) 


* 


Since  x  is  we  are  explicitly  assigning  equal  weights  to  the  spin 

rate  and  coning  rate.  As  the  moment  is  a  response  to  the  coning  motion,  the 

appearance  of  is  natural.  The  primary  characteristic  of  the  high  Reynolds 

number  Stewartson-Wedemeyer  theory  is  the  existence  of  a  dominant  inertia  wave 
mode.  This  resonance  mode  is  caused  by  an  interaction  of  the  spinning  liquid 
with  the  coning  motion  and,  thus,  the  spin  rate  is  as  important  as  the  coning 
rate.  At  very  low  Reynolds  numbers,  these  inertia  waves  are  strongly  damped 
and  the  influence  of  spin  on  the  liquid  moment  should  be  much  weaker.  Indeed, 
Miller8  found  no  dependence  of  the  roll  moment  on  spin  for  t  <  .3  and  very  low 
Reynolds  numbers.  If  the  roll  moment  is  independent  of  spin,  the  liquid  roll 
moment  coefficient  is  proportional  to  t  (C^r^  =  x  CR). 


The  frozen  liquid  value  (Re  =  0)  of  the  liquid  in-plane  moment  coeffi¬ 
cient  for  constant  amplitude  coning  motion  is10 


CLIM  =  1/2  -  x  [1/4  +  c2/3a2] 


(2.7) 


Thus  the  in-plane  moment  must  always  be  strongly  dependent  on  spin.  The  side- 
moment  coefficient,  however,  for  t  <  .3  and  low  Reynolds  number  is  probably 
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proportional  to  t  (C^m  =  t  C$)  and  for  these  conditions  both  and  CR  may 

be  more  appropriate  non-dimensional  forms  of  these  moments.  (These  new 
coefficients  may,  of  course,  be  functions  of  t.) 

III.  ROLL/SIDE  MOMENT  RELATION 
In  Reference  5,  it  is  shown  that 


or 


cLrm  =  -CLSM  +  (™/2)  E1  -  (4/3)  (c/a)2^ 

CR  =  -Cs  +  (e/2)  [1  -  (4/3)  (c/a)2] 


(3.1) 

(3.2) 


These  equations  are  valid  for  linear  fluid  mechanics  and  state  that  for  the 
simple  case  of  constant  amplitude  motion  (e  =  0)  the  liquid  roll  moment  is  the 
negative  of  the  liquid  side  moment.  Since  no  boundary  layer  assumptions  were 
used  to  obtain  Eqs.  (3. 1-3. 2),  they  should  apply  to  the  numerical  results  of 
Reference  12.  If  the  complete  nonlinear  fluid  equations  are  used,  additional 
significant  contributions  to  the  quadratic  liquid  roll  moment  may  be  present. 
It  is  interesting  to  note  that  Eq.  (3.1),  for  constant  amplitude  coning 
motion,  is  the  linear  condition  for  the  liquid  moment  about  the  trajectory  to 
be  zero. 


MLX  =  MLX  cos 


“t 


MLSM  sin  “t  =  MLX  + 


MLSM  K1 


(3.3) 


Thus  the  ratio  of  the  magnitudes  of  the  roll  moment  to  the  side  moment  is 
tan  ctj. . 


Recently  Pope13  very  carefully  analyzed  the  yawsonde  records  of  twelve 
unstable  liquid-filled  projectiles.  All  projectiles  showed  a  growing  coning 

motion  and  a  rapid  despin.  Two  of  these  missiles  carried  low  viscosity 
liquids  (Re  ~  106  ),  four  carried  silicon  oil  (Re  ~  40),  five  carried  corn 
syrup  (Re  ~  10),  and  one  had  felt  wedges  impregnated  with  liquid  white 
phosphorus. 

Since  the  aerodynamic  roll  damping  is  easily  measured.  Pope  was  able 
to  obtain  liquid  roll  moment  coefficients.  He  was  unable  to  separate  the 
aerodynamic  and  liquid  side  moments  in  the  observed  motion  and,  therefore,  ob¬ 
tained  "equivalent"  side  moment  coefficients  (C^eq)  containing  both  contribu¬ 
tions.  He  was  able  to  make  estimates  of  the  aerodynamic  side  moment  coeffi¬ 
cients  (CSMA).  These  three  coefficients  are  given  in  Table  1  for  the  twelve 

yawsonde  records.  The  rough  agreement  of  C|_RM  and  C$M£q  with  Eq.  (3.1)  is 
very  gratifying.  w 

Finally,  Miller's  liquid  roll  moment  data6  for  spinning  cylinders  con¬ 
ing  at  a  constant  cone  angle  of  20°  are  plotted  in  Figure  1  as  open  circles. 
Clsm's  are  computed  by  the  linear  boundary  layer  theory  of  Reference  10  and 
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plotted  as  solid  circles.  The  qualitative  agreement  for  the  boundary  layer 
theory  at  such  low  Reynolds  numbers  is  quite  encouraging. 


Table  1.  Pope's  Results 


Yawsonde 

Re 

T 

CLRM 

CSMEQ 

CSMA 

404  « 

2  x  106 

0.090 

-0.05 

0.05 

<0.01 

1339 

1.8  x  106 

0.090 

-0.02 

0.04 

-0.005 

1866 

45.2 

0.123 

-0.055 

0.050 

-0.004 

1867 

45.2 

0.123 

-0.055 

0.050 

-0.004 

1868 

45.2 

0.123 

-0.060 

0.053 

-0.004 

1955 

20 

0.087 

-0.04 

0.04 

-0.010 

1293 

10 

0.090 

-0.03 

0.02 

-0.008 

1313 

10 

0.088 

-0.02 

0.025 

-0.008 

1585 

10 

0.091 

-0.025 

0.02 

-0.008 

1587 

10 

0.093 

-0.04 

0.02 

-0.008 

1588 

10 

0.095 

-0.03 

0.02 

-0.008 

1693  b 

0.094 

-0.025 

0.025 

-0.005 

a  Fill  ratio  for  404  was  87%,  all  others  were  100%. 

b  This  projectile  contained  white-phosphorus-impregnated  felt  wedges.14 


14.  D'Amico ,  W.  P. ,  " Aeroballistic  Testing  of  the  XM825  Projectile:  _  Phase 
III ,  High  Muzzle  Velocity  and  High  Quadrant  Elevation US  Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground ,  Maryland ,  BRL  Memorandum 
Report  ARBRL-MR-03196 ,  September  1982.  (AD  B068511) 
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Figure  1.  CLRM  From  Gyroscope  Tests  (Ref.  6),  c/a  =  4.291. 
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